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Трёхмерная контактная задача для трансверсально изотропного тела” 
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Изучена пространственная контактная задача с неизвестной областью контакта для трансверсально- 
изотропного упругого полупространства, граница которого перпендикулярна плоскостям изотропии. ЖЕСТ- 
кость границы упругого тела зависит от направления, поэтому для кругового в плане штампа область кон- 
такта, как правило, не является круговой. Задача сведена к интегральному уравнению (ИУ) относительно 
контактного давления, ядро которого не содержит квадратур. Для решения ИУ применяется численный ме- 
тод Галанова, позволяющий одновременно определить область контакта и давление в этой области. Простой 
вид ядра ИУ облегчает его регуляризацию, параметр которой зависит от шагов сетки и параметров анизо- 
тропии. Для отладки программы использовано известное точное решение задачи для штампа в форме эл- 
липтического параболоида. Сделаны расчёты для разных трансверсально-изотропных материалов при внед- 
рении конического и пирамидального штампов. 

Ключевые слова: теория упругости, контактная задача, трансверсально-изотропное полупространство, 
метод Галанова. 


Введение. Механические свойства трансверсально-изотропных материалов вызывают интерес у 
исследователей, поскольку касаются ряда важных материалов, имеющих гексагональную кри- 
сталлическую решётку [1]. Развиваемый метод позволяет эффективно оценивать твёрдость и 
контактную прочность материалов, приповерхностные свойства которых могут существенно зави- 
сеть от направления. Ядро интегрального уравнения контактной задачи ранее было получено в 
виде двойного интегрального преобразования Фурье [2]. Затем при помощи теории обобщённых 
функций удалось представить это ядро в виде свободном от квадратур [3]. Такой вид ядра, ввиду 
простоты его регуляризации в особых точках, сделал возможным применить для решения кон- 
тактной задачи с неизвестной областью контакта метод нелинейных граничных ИУ типа Гаммер- 
штейна, развитый Галановым [4, 5]. Твёрдость по Бринеллю для эллиптического штампа может 
быть оценена на основе точного решения для эллиптического параболоида [2]. В представленной 
работе сделаны расчёты для разных материалов при внедрении конического штампа и штампа в 
форме четырёхугольной пирамиды (твёрдость по Виккерсу). Определены области контакта, дав- 
ления и значения вдавливающей силы при заданной осадке штампа. 

Постановка задачи. В декартовых координатах рассмотрим трансверсально-изотропное упругое 
полупространство х> 0, граница которого перпендикулярна плоскостям изотропии 2 = сопв(. За- 
кон Гука, включающий пять упругих параметров, и уравнения равновесия приведены в [2, 3]. 
Пусть при х = 0 в полупространство внедряется абсолютно жёсткий штамп, основание которого в 
области контакта описывается функцией Ку,2). Штамп вдавливается без перекоса центрально 
приложенной силой Р, испытывая осадку б. При заданных упругих параметрах, значении дб и 
функции Ку,2) требуется определить область контакта ©, контактное давление 4(у,2) и силу Р. 
На основании решения задачи Буссинеска [1] ИУ контактной задачи можно записать в форме 


Да(’,>)К(у-У»2-2)ву ав, =56-Г(у,?), (у,2)е ©, (1) 


где ядро представляется двойным интегралом Фурье (3.4) [2] или (1.4) [3]. 





* Работа выполнена при финансовой поддержке РФФИ (грант 12-01-00065). 
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При помощи теории обобщённых функций ядро удаётся представить в форме свободной 
квадратуры [3]: 
т -т,)у? у? 
(т) У м РС 
2пАз р 
р = т,В25, - т,В2б, -4 (т, - т,) 22. 65., (2) 
Ану С А а 2 А 
— „В = +1 2 [=1,2), р Е 
м А, 3+ Аз ' (т }\зу ме } ь А. 
Здесь Аи, Алз, Азз, Ада, Ав — упругие параметры, у, у> удовлетворяют уравнению 
У‘А Ал ыы у’ [А 1153 и А,з(А 13+ 2А ) Аз Ак =0. (3) 
В табл. 1 (колонки 3—7) приведены параметры анизотропии (2), (3) для ряда материалов, пара- 
метры упругости которых А’ табулированы в [1]. 
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Таблица 1 
Параметры анизотропии (безразмерные) 

№ Материал У у Уз т, т, Ч И 

1 | А>Оз 2,786 0,3974 0,8890 4,543 0,2201 0,7283 0,7248 
2 |0 3,269 0,3568 1,061 5,206 0,1921 0,5701 0,6504 
3 Мд 2,050 0,5041 0,9791 2,782 0,3595 0,6778 0,6950 
4 |9С 2,859 0,3808 0,7786 5,994 0,1668 0,8911 0,8304 
5 п 1,759 0,6324 1,327 2,066 0,4840 0,5204 0,6384 
6 | (45 3,186 0,3432 0,9198 4,219 0,2370 0,5849 0,6167 
7 | Углеволокно 7,492 1,568 4,790 4,132 0,2420 0,2941 0,6671 
8 | Графит 105,1 0,0003322 | 0,0007955 7256 0,0001378 22,95 0,8548 
9 | Сапфир 2,336 0,4293 0,8848 3,921 0,2550 0,7695 0,7523 
10 | Древесина (ель Дугласа) 13,79 0,1227 0,7101 76,22 0,01312 0,5306 0,6021 
11 | Керамика Р7Т-4 1,198 0,6907 0,8393 1,415 0,7069 0,7066 0,6406 
12 | Композит (60 % волокон) 22,32 0,4745 1,821 25,91 0,03859 0,4099 0,6602 
13 | Бедренная кость человека 4,157 0,2787 0,9429 4,443 0,2251 0,4779 0,5412 
































Поверхность полупространства, в зависимости от материала из которого оно состоит, мо- 
жет быть более жёсткой как в направлении оси у, так и в направлении оси 2. При действии на 
границе полупространства в начале координат нормальной сосредоточенной силы Р, для нор- 
мальных перемещений на осях уи 7х можно получить формулы [3] 








Р, и, (т, т, уу, 
и, (0,у,0) = ——.—^, и, = , 
ле и У} | т, (т, +1} у, т, (т, +1), | 
я И (т, — т, Уз 


и, (0,0;2)= 











х 
214, 2’ “= 2| (т, т, У — ту? + т»\! | | 

Значения безразмерных величин Шу Ц. приведены в двух последних колонках табл. 1 и характе- 
ризуют нормальные перемещения точек поверхности, лежащих на координатных осях и равно- 
удалённых от начала координат, где действует сосредоточенная сила. Видно, что для материалов 
1, 4, 8, Зи 11 поверхность полупространства жёстче в направлении оси 2 (соответствующее пе- 
ремещение меньше). Для остальных материалов поверхность жёстче в направлении оси у. 


Направления экстремальной жёсткости границы полупространства после перехода к по- 
лярным координатам р, а можно определить из формул [3] 








и, (0, рсоза, рэта) _ т, и (а) — (т, = т, ) Уз в С: (а), (а) р 
р (а) = т,А; (а)с, (а)- т»? (а)с, (а) - 4 (т, - т, ) п? аб, (а), (а)с, (а), (4) 
В, (а) = (т, +1)уз со5° а+2 т? а, С,(а)= № 05? а + п? а. 
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В табл. 2 для материалов из табл. 1 даны значения углов а а при которых функция и(а) 


тт / тах / 


(4) достигает на отрезке [0; п/2] соответственно наименьшего и наибольшего значения. Эти углы 
характеризуют направления наибольшей и наименьшей жёсткости поверхности полупространства. 














Таблица 2 
Направления экстремальной жёсткости 
Материал 1 2 3 4 5 6 7 8 9 10 11 12 13 
О и п/2 0 0 п/2 0 0 0 п/2 п/2 0 п/2 0 0 
Ч тах 0,699 п/2 1,20 0 п/2 1,08 п/2 0 0,503 | 0,833 0 п/2 1,21 












































Как видно, экстремальные направления безразмерного перемещения (4) часто (но не все- 
гда) совпадают с направлениями осей координат у, 2. 
Решение задачи. Для решения ИУ (1) при условии 9(у,2) = 0, (/,2)>Е©, используем метод нели- 
нейных граничных ИУ [4, 5], позволяющий одновременно определить область контакта, контакт- 
ное давление, нормальное перемещение материала вне области контакта. Предположим, что об- 
ласть контакта целиком содержится в прямоугольнике $ = [и <В,, |2| < а}. Прямоугольник 


$ покроем равномерной сеткой из л7 узлов с шагами /, по оси уи №» по оси 2. При расчёте значе- 

ний ядра в этих узлах его особенности сглаживались по формулам 

ть (5) 
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Регуляризация (5) обеспечивает сходимость метода и отладку программы, давая хорошее совпа- 

дение с точным решением для эллиптического параболоида [2], когда 


У(У- У.) +п(2-25} > У(У-У) +1(2-2,} + (\+п5., 5. = 


о - 
Ре 28 "28 (6) 


Для случая (6) решение ИУ (1) с ядром в форме двойного интеграла Фурье получено в виде [2] 


2 2 
2 
9(у,2) = у = Р= [а(у,2)вуа = 5`аь» (7) 
о 


(т, — т, )уз 226, (9) 5, (9) соз? 94 














б=а4 ‚ (8) 
. 8А; ] р(6)г (9) 

Ь? а? в К _ С _ а . 2 А 

5. +. 5, ка’ г (0) {2 С05? 9 + $? 0, (9) 


где использованы обозначения (4) и 
те 215, (9)5, (9) со5* 04 ре |< (9), (9) соз? 0 п? ва@ 
Гео ‘Г 9) 

При заданных величинах 5, А!, К, отношение полуосей эллипса контакта а/Б определяется из вто- 
рого соотношения (9). Затем величина а находится из первой формулы (9), величина 4 — из (8). 
Вдавливающая сила рассчитывается по второй формуле (7). 

При расчётах брали а = № (прямоугольник $ — квадрат), случаи материалов соответ- 
ствуют табл. 1. Число узлов бралось от 81 до 289. 

Использовались безразмерные обозначения (штрихи далее опускаем) 














К Г 
ры а бе. К, = т, К, = 2, а И. 
4 а 4 4 а а а 
, , , Ч(Уу,2 ' 9 , Р 
К = , 9'(у’,2') = } ЧФ = °, Р 2 
Аз Аз Азваь 
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Отметим, что при вдавливании эллиптического параболоида может возникать круговая область 
контакта (К = 1). При вдавливании кругового параболоида (А, = А), как правило, возникает эл- 
липтическая область контакта. При этом оказывается, что для материалов 1, 4, 8, 9 и 11, когда 
поверхность полупространства жёстче в направлении оси 2, область контакта вытянута вдоль 
ОСИ 2. Для материалов 2, 3, 5—7, 10, 12 и 13, когда поверхность тела жёстче в направлении оси у, 
зона контакта вытянута вдоль оси у. 

Аналогичный вывод можно сделать для штампа в форме кругового конуса, когда 


Е(у,2) = У? +2°, 
и для штампа в форме правильной четырёхугольной пирамиды, когда 
Г(у,2) = тах (|у|,[2|). (10) 
На рис. 1 схематично показаны узлы сетки, принадлежащие области контакта, при вдавливании 
пирамиды (10). При этом брали сетку из 13х13 узлов, д = 1. Как видно, по сравнению с изотроп- 
ным материалом, зона контакта увеличивается в том направлении, в котором поверхность тела 


более жёсткая. Из материалов, приведённых в табл. 1, наименьшая площадь области контакта 
наблюдается для графита. 






































































































































































































































пинииниини о. са св 
онииинининио ив ни ив 
ини инаи ив Г ив 
нии ини пинниошншна 
нии инаи ив Г Г 
ониинииннио ив ии ив 
ини наи пси пси 
а) 6) 


Рис. 1. Узлы сетки в области контакта: а — к затемнённым квадратикам для изотропного материала добавляются 
незатемнённые для материалов 4, 11, ромбики для материалов 2, 5, 10 или ромбики и кружки для материала 12; 
6 — незатемнённые квадратики для графита, все квадратики для углеволокна. 


Заключение. Решена новая пространственная контактная задача с неизвестной областью кон- 
такта для трансверсально-изотропного тела. В отличие от случая, когда плоскости изотропии па- 
раллельны границе полупространства [1], здесь для вытянутого в определённом направлении 
штампа форма области контакта существенно зависит от ориентации штампа. 
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ТНВЕЕ-ОТМЕМ$ТОМАЕ СОМТАСТ РВОВЕЕМ ЕОВ А ТКАМ$УЕВ$Е1У 
ТбОТВОРТС $011” 


О. А. РотНаг$Юу, О. В. Ваубуап 
(Роп Зае Тесптка! Упмег$Ку) 


Тре рава! сощасЕ ргоМет ийёН ап ипкпоигп сог(асЕ аота! 1/5 пуезВдейеа ГГ а Егапзуегсе/у 1боторгс е/аз с Ва/!Е- 
5расе те Боипаагу оЁ иРЕСВ 15 регрепаси/аг Ю те р/апез о! 15оёгору. Рог а сгси/аг рипсй, пе согвасеЕ гопе, а а 
ИЕ, 15 ПОЁ а стае Бесаизе {те 5Итез5 оЁ те е/аз8с зо Боипаагу дерепаб оп йе а!гесйоп. Тпе ргоМет (5 ге- 
исеа © ап И\едга! едиавоп (ТЕ) и" гезресЕ 0 {те сощасЕ ргезиге #е Кетае! о! и/сВ аое$ поЕ тс/иаде диадга- 
{игез. байапоу5 питейса/ теоа иИсВ такез # роз {0 а ептте хтийапеоиз/у пе сопЁасЕ гопе апа пе 
сог(асЕ ргеззиге 15 изеа ю 5о№е {те ГЕ. Тре ятр/е Югт о! {пе [Е Кете! аси гедшатета й Бу изтд а рагате- 
{ег ир/сй Черепа$ оп тезВ И\егуа/5 аз ие! аз оп атб5огору рагатеЁегз. А ие/-кпоилт ехасЕ зо/иНоп Ю а рипсй 
т {те Югт оЁ ап еМрёса! рагабо/о/А 15 изей юЮ ету {те сотршег ргодгат. ТВе питепса! апа/ух5 Ваз Бееп 
таае ог ЙТегепЕ {апзуегсе/у 5оторс теепа5 сощасйпд ит согиса! апа ругат!аа! рипсйес. 

Кеуигогаб: {Пеогу оЁевзйсйу, сог(асЕ рго ет, гапжуегсе/ Г5оторгс Ва!5расе, ба/апоу5 те{тоа. 
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